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We propose the measurement of the elliptic flow of hadron resonances at the Relativistic Heavy 
Ion Collider as a tool to probe the amount of hadronic final state interactions for resonances at 
intermediate and large transverse momenta. This can be achieved by looking at systematic devi- 
ations of the measured flow coefficient V2 from the scaling law given by the quark recombination 
formalism. Our method can be generalized to explore the structure of exotic particles, such as the 
recently found pentaquark O + (1540). 
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Recently, progress has been made in our understand- 
ing of the hadronization of a quark gluon plasma created 
in high energy heavy ion collisions. Both the yields and 
the anisotropic flow v 2 of hadrons at transverse momenta 
above 2 GeV, observed at the Relativistic Heavy Ion Col- 
lider (RHIC), are described well by a combination of two 
mechanisms: (i) the statistical recombination of thermal 
constituent quarks from a quark gluon plasma phase dur- 
ing the phase transition, (ii) fragmentation of jets which 
experienced energy loss due to the surrounding medium 

The elliptic flow v 2 is defined as the second Fourier 
coefficient of the azimuthal hadron spectrum [(| 



d 2 N 



dN 



[l + 2v 2 {P T ) cos 2$ + ...]. (1) 
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It measures the elliptic anisotropy in the spectrum since 



v 2 (P T ) = (cos 2$) 
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Such an asymmetry naturally arises in heavy ion colli- 
sions with non- vanishing impact parameter b > 0. 

It has been shown that in the region where recombina- 
tion of partons dominates the hadronization process — 
i.e., for Pt < 4 (6) GeV for mesons (baryons) — and 
mass effects are suppressed — i.e., Pt > 2 GeV — v 2 
obeys the simple scaling law, 



v 2 (Pt) = nv^Px/n) 



(3) 



where n is the number of valence quarks of the hadron 
0- @- We will assume throughout the paper that 



strange quarks and light quarks have the same elliptic 
flow v\. Measurements of strange hadrons support this 
assumption 2]. At higher transverse momentum, Pt > 6 
GeV, where fragmentation of partons dominates, we ex- 
pect the v 2 of all hadrons to lie on a universal curve since 
the effect of fragmentation functions largely cancels in 
Eq. (0) . We refer the reader to Refs. 0, for a more 
detailed discussion. 

The scaling law © has been impressively confirmed by 
measurements at RHIC. Pions and kaons (n — 2), pro- 
tons, A and S (n = 3) all fall on one universal curve, if 
V2 and Pt are divided by the number n of valence quarks 
[tj, El • Slight deviations can be observed and are largest 
for pions. This can be due to its nature as a Goldstone 
boson — its mass is much smaller than the sum of its 
constituent quark masses — or because a large fraction 
of all pions are not created at hadronization but by the 
subsequent decay of hadron resonances 0] . This poses 
the interesting question how the existence of a hadronic 
phase generally affects hadron production at intermedi- 
ate and large momenta. We suggest that the amount of 
hadronic rescattering can be determined by measuring v 2 
for several hadronic resonances. 

Hadrons can experience rescattering after hadroniza- 
tion, until kinetic freeze-out takes place. Experimentally, 
it has been found that the kinetic freeze-out temperature 
at RHIC ^| i s much lower than the expected QCD phase 
transition temperature |13|. Resonances produced from 
a hadronizing quark gluon plasma decay in the hadronic 
medium. If their decay products rescatter, the signal is 
lost in the experimental measurement, since the origi- 
nal resonance cannot be reconstructed from a correlation 
analysis of its daughter hadrons. Likewise, hadrons from 
the medium may scatter into a resonance state and thus 
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contribute to the final measured yield. The total mea- 
sured resonance yield will therefore have two contribu- 
tions: 

1 . primordial resonances: resonances produced from a 
hadronizing quark gluon plasma whose decay prod- 
ucts have not rescattered (QGP mechanism). 

2. secondary resonances: resonances produced in the 
hadronic final state via hadron-hadron rescattering 
(HG mechanism). 

In this paper, we will show that vi differs between the 
two cases. Key to our analysis is the observation that the 
elliptic flow i>2 will be additive for any type of composite 
particle with respect to the v 2 of its constituents. Note 
that this feature is mimicked in a hydrodynamic model 
when mass effects are neglected. Then the elliptic flow 
of all particles follows the same line through the origin, 
V2{Pt) oc Pt- The scaling law simply maps this line into 
itself. Only the deviation from the ideal hydro behav- 
ior with increasing Pt, manifest in the saturation of i>2 
at intermediate Pt, makes the scaling law non-trivial. 
The V2 of a Kq meson in the recombination domain will 
therefore be the sum of the v 2 contributions of the d and 
s quarks if it has been formed from a hadronizing QGP, 
or the sum of the contributions of the K + and 7r _ if it 
has been formed through the coalescence of a kaon and 
a pion in the hadron phase. Our discussion of deriv- 
ing the formalism for the coalescence of partons can be 
extended to the coalescence of hadrons in a straightfor- 
ward way 111 LI La] Therefore we can here assume the 
same additive feature for hadron-hadron coalescence as 
well. Naively, if the kaon and pion themselves are formed 
by quark recombination, we expect a scaling of the Kq 
with n = 4 in the HG case, while it scales with n = 2, 
like the stable mesons, if it is created at the phase tran- 
sition. The scaling will be altered by the possibility that 
one or both of the coalescing hadrons could come from 
fragmentation, which breaks the scaling law. The final 
result will be a mixture of both scenarios. 

We assume that the V2 of stable hadrons (stable with 
regard to strong interactions) is not affected by hadronic 
interactions for Pt above 2 GeV. This is a sound as- 
sumption in light of the experimental facts, as discussed 
above. Let r(Pr) be the fraction of resonances that have 
escaped from hadronization without rescattering of the 
decay products vs. the total measured yield. Given the 
knowledge of V2 in the two limiting cases, the QGP mech- 
anism and the HG mechanism, the fraction r(Pr) can be 
determined by measuring V2 at intermediate Pt , between 
2 and 6 GeV via 

^measured = ^p^QGP + ^ _ r (p T ) )l; HG^ (4) 

One may be able to deduce the widths and cross sec- 
tions of resonances in the hadronic medium by measuring 
r(Pr) for several resonances. Our scheme can be gener- 
alized to explore the structure of exotic particles, e.g. a 



large molecular state, as discussed in the case of the re- 
cently discovered pentaquark ©+(1540). We will discuss 
this in more detail below. 
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FIG. 1: (Color online) Elliptic flow V2 as a function of Pt 
for the £*(1385) (top) and the K<,(892) (bottom). Dashed 
line: resonance contribution from hadronization (QGP); dot- 
ted line: resonance contribution from coalescence in a hadron 
gas (HG); solid line: calculation for a 30% contribution of 
hadronic resonance formation. 

Figure ^ shows the elliptic flow v 2 as a function of 
transverse momentum Pt for the S*(1385) (top) and the 
Kg (bottom) resonances. The dashed lines refer to pri- 
mordial resonance production from a hadronizing QGP 
via parton recombination at low and intermediate Pt 
and fragmentation at large Pt- The dotted lines refer 
to these resonances being produced via hadronic final 
state interaction through hadron-hadron recombination: 
A + 7r — > S* and K + + ir~ — > Kq. These two scenarios 
establish lower and upper limits for V2- The solid line 
represents a calculation with contributions from both, 
the QGP and the HG mechanisms. For demonstration 
purposes we have here chosen r(-Pr) to be constant, with 
a value of 70%. r(Pr) is the quantity that can ultimately 
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be extracted from data. 

For quark recombination we used the scaling law (0) 
together with the parton v\ extracted and discussed in [3, 
3j. Since the fragmentation functions for the resonances 
are unknown, we utilize the universal curve for the vi 
contribution from fragmentation as obtained in 0, . 

Due to the transition from the recombination to the 
fragmentation regime, the maxima of the HG and QGP 
curves do not scale like 5/3 and 4/2 for the S* and K^ 
respectively. Instead, the HG curve is brought down by 
contributions from pions, kaons or As from fragmenta- 
tion, violating the scaling law. Still, the difference be- 
tween the two limiting cases can be as large as 30% be- 
tween 2 and 6 GeV. 

The fraction of the total resonance yield being pro- 
duced by re-generation in the hadronic phase cannot be 
inferred through the final measured resonance yield since 
the amount of signal loss in resonance reconstruction due 
to rescattering of the resonance decay products cannot be 
experimentally determined. The scheme outlined above 
to determine r{Pr) therefore provides complimentary in- 
formation to the resonance yield measurements via in- 
variant mass reconstruction. The determination of r(Pr) 
allows for an estimate of the gain term of resonance pro- 
duction in the hadronic phase - together with the yield 
information this allows for a far better estimate of the ac- 
tual QGP production of resonances and for an improved 
estimate of the signal loss due to rescattering of the de- 
cay products. These estimates are of high relevance for 
the interpretation of statistical model fits with respect to 
resonances yields. 

We can now utilize the scaling law for the elliptic flow 
- and its transition to a universal value for V2 in the 
fragmentation domain - to probe the structure of exotic 
particles, such as the 8 + (1540) pentaquark state. The 
discovery of this novel hadronic state with five valence 
quarks (uudds) has recently been reported in 0|. Its 
existence was predicted in 1997 by Diakonov et al. in the 
chiral quark soliton model (CQSM) |l5j. Several other 
gro ups have by now confirmed this important discovery 
[l6L fl7l fl8| . The 6 + is the first unambiguous five quark 
state, since other candidates for five quark states, such 
as the A(1405), have quark compositions that can mix 
with three quark states like uddds <-> uds. Recently, 
NA49 reported the observation of two more states that 
could be the two missing pentaquarks in the antidecuplet 
[T^ | . The yield of the 8 + in heavy ion collisions has been 
estimated in p(ll2lT|. 

While the discovery of the + has undoubtedly es- 
tablished the existence of pentaquarks, its structure, and 
even its spin and parity, have not yet been experimentally 
verified. It is still an open question whether pentaquarks 
like the + have a compact structure with all five valence 
quarks being closely confined, just as the three quarks in 
the nucleon, or whether they are molecular bound states 
of a baryon and a meson. A bound state of two diquarks 
and an antiquark was recently suggested as well 22] . 

Loosely bound states with large geometric cross section 



whose binding energies are much less than the temper- 
ature, like the deuteron, are dissociated in the hadron 
phase even if they were produced at hadronization. Such 
a medium induced breakup reaction does not allow for 
the reconstruction of the state via two-particle mass spec- 
tra, since the initial momentum of the scattering partner 
is unknown. These bound states are unobservable. Thus, 
the loosely bound states and molecular states observed 
in heavy ion collisions are produced close to the kinetic 
freeze-out hypersurface |2l| , in the case of the pentaquark 
by K + +n — > + and K° +p — » + (summarily denoted 
as KN). The hadronic phase can be viewed as a filter 
that tends to let the + from hadronization only survive 
if it is a compact five quark state. 

With respect to the azimuthal anisotropy of pen- 
taquarks produced in heavy- ion collisions at RHIC, we 
expect it to fulfill the scaling law © with n = 5. This 
would be an impressive manifestation of its exotic five 
quark character. On the other hand, we expect modi- 
fications to hold, as we discussed them for conventional 
resonances. However, we want to point out, that in the 
case of conventional baryon resonances the scaling law 
holds for n — 3 in the QGP scenario, while it is roughly 
n = 5 in the HG case. For a pentaquark it would be 
n = 5 in both cases. The difference would merely come 
from the possibility that for the reaction N + K — > + 
one or both of the initial hadrons could come from frag- 
mentation. The same is true for a deuteron, that scales 
approximately with n = 6 in both cases. 

If the + is a molecular state, the momentum fraction 
xk of the kaon in a lightcone frame is expected to be 
approximately given by tk = Mk/(Mk + Mjy) ~ 0.35, 
where and Mjv are the masses of the kaon and nu- 
cleon, respectively. We therefore use the wave functions 

\4>q(xk,xn)\ 2 = S(x K -0.35), (5) 
\cj>d{x pi x n )\ 2 = S(x p - 0.5). (6) 

In Fig. [21 we show vi for the different scenarios dis- 
cussed above. Note that the scenario with all pen- 
taquarks from hadronization (QGP) is now providing the 
upper limit, while it was the lower limit in the case of 
conventional resonances. We include the deuteron in the 
study to give an example of a well known molecular state. 
Above Pt — 3 GeV we find a considerable difference, up 
to 20%, between the scenario where the pentaquark is 
created at the phase transition and is therefore a com- 
pact five quark state — and the scenario where it origi- 
nates from KN coalescence. If a KN coalescence takes 
place to form a + between 4 and 8 GeV, it is likely that 
either the kaon or neutron (or both) stems from fragmen- 
tation which lowers the resulting V2- At very high Pt, 
both the kaon and the nucleon have to come from frag- 
mentation. This results in a power law spectrum that 
falls off like P^f 2b , approximately twice as rapidly as a 
directly fragmenting particle with ~ P^ b . Therefore, di- 
rect fragmentation will dominate at very large Pt and 
will eventually lead v-i down to the universal curve. 
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FIG. 2: (Color online) Pt dependence of V2 for the pen- 
taquark + (top) and the deuteron d. Thick solid line: frag- 
mentation and recombination of a 5 (6) quark state for + 
(d); thick dashed line: fragmentation and recombination of 
K + N (p + n) for Q + (d) as true for a molecular state; thin 
solid line: direct fragmentation (universal curve) ; thin dashed 
line: recombination of K + N (p + n) only. 



For the deuteron as a shallow molecular bound state we 
expect pn recombination at the kinetic freeze-out to dom- 
inate. The deuteron is therefore ideal as a benchmark 
measurement. First results for the v-i of the deuteron 
have been announced by the PHENIX collaborations . 
As expected, this particle seems to follow the scaling law 



with n = 6 within the large error bars. With the upcom- 
ing high statistics Au+Au run, error bars and Pt cover- 
age should be greatly improved, allowing for a gauging 
of the formalism. 

Let us comment on the third possible scenario for the 
+ , that was advocated by Jaffe and Wilczek j2^|. If it is 
a collection of two tightly bound (ud) diquarks and an s 
quark, recombination would predict the Pt dependence 
of «2 to be very similar to that of a "democratic" 5 quark 
state. This, however, is only true if the diquarks are 
also formed at hadronization and the individual quarks 
have not been correlated early on with enough time to 
have thermalized diquark states. In that case one would 
expect a scaling with n = 3. 

In summary, we propose the measurement of the ellip- 
tic flow of hadronic resonances at the Relativistic Heavy 
Ion Collider. After recombination scaling has been con- 
firmed for the vi of stable particles, deviations from this 
scaling for resonances can directly measure the fraction 
of resonances that escape from hadronization without 
rescattering in the hadronic phase at intermediate and 
large transverse momentum. This allows conclusions to 
be drawn about the importance of the hadronic phase in 
this region of phase space. 

In turn, this method can be applied to explore the 
structure of exotic particles, such as the recently found 
O + (1540). The amount of rescattering in the hadronic 
phase is sensitive to the size of the system, eventually 
distinguishing between a compact hadron and a molecu- 
lar state. The deuteron, with a known molecular struc- 
ture, can serve as a precise benchmark for these measure- 
ments. The same method is also applicable to determine 
the structure of other exotic candidates with a possible 
molecular structure such as A(1405), ao(980), /o(980), 
and so on gHUHl. 
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